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ABSTRACT 

Electromagnetic decays of the scalar mesons are shown to be constrained by chiral 
^< ■ symmetry as a consequence of the fact that, in the chiral limit, the two and three- 

Q\ \ point functions < SS—PP > and < WS > satisfy super-convergent dispersion 

CNI ■ relations. The QCD asymptotic behavior of the latter is canonical and it can be 

saturated by a finite number of resonances. The corresponding chiral lagrangian 
for vector and scalar resonances is constructed. Matching to the correct asymptotic 
| structure generates non-minimal terms which have so far been ignored. It is found 

I/-) ■ that the width a (980) — > 27 can be naturally reproduced suggesting that the 

CO | do (980) is not an exotic particle. 

■"xl" ! 

1- Introduction: 

The chiral symmetry limit of QCD, which is expected to be a reasonably good 
r£3 ' approximation in the light quark sector, may be at the origin of some of the surprising 
( at first sight) properties that the scalar mesons seem to have. We will concentrate 
Ph. here on the question of the two photon width of the a (980) (since this is a gamma- 
gamma meeting) and it will be shown that the fact that this width is anomalously 
small (as compared, for instance, with the expectation from the naive quark model) 
has a natural explanation in this limit, and is by no means incompatible with a 
standard nonet status for the do- 

Our starting point is to assume exact chiral symmetry (i.e. we set m u = 
m.<2 = Tn s = in the lagrangian) and , furthermore, we will work in the leading 
N c approximation. In these limits, QCD still cannot be solved in closed form, but 
analytic approximation methods, at least, are available. The one that we will use 
is similar to the boundary- layer technique used, for instance, in hydrodynamics and 
consists in matching three asymptotic expansions for a given correlation function: a) 
For large momenta, the operator product expansion (OPE) generates an expansion 
in inverse powers of the momenta b)for small or moderate momenta the correlation 
function is dominated by the contributions from the light resonances and c)at very 
low momenta the expansion is generated by chiral perturbation theory (xpt ). The 
method is effective only for those n-point functions which transform in a non trivial 
way under the chiral group and thereby satisfy super- convergent dispersion relations. 
Only in that case is it possible to match the QCD asymptotic behavior with a finite 
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number of resonances. This method has been used in the past for a diversity of 
applications, for the first time by Weinberg who obtained an estimate for the mass 
of the ai resonance. Das et al0 have evaluated in this way the 7r + — 7r° electromagnetic 
mass splitting. More recently, similar ideas where discussed in the context of the 
effective lagrangian for resonances 0] and an interesting application to the question 
of whether the mass of the u quark can be zero can be found in ref . . 
2- Meson spectroscopy: 

In the exact chiral limit, the spectrum of the lightest resonances (mesons as 
well as baryons) in the QCD spectrum is expected to consist of exactly degenerate 
octets (which, for the mesons, become nonets in the leading N c approximation). This 
is not what is observed in practice, of course, but already first order perturbation 
theory in terms of the quark mass hamiltonian predicts a pattern of splitting, ideal 
mixing, Gell-Mann-Okubo mass formulas, which is satisfied, in practice, to a fairly 
good approximation. The fact that the pion nonet does not obey ideal mixing is due 
to the axial anomaly and is well understood in terms of the large N c expansion and 
Xpt techniques (e.g. f| ). The only major exception to these rules seem to be the 
scalar resonances. The present candidates to form a nonet that one can find in the 
review of particle properties |J are: 

1 = 1: a (980), 1 = 1/2: #£(1430) 

1 = 0: / (975), 1 = 0: /„(1400) (1) 

The decay properties of the / (975) and the / (1400) strongly suggest that the flavor 
transformation properties of the former is essentially like ss and the latter is uu + 
dd. The pattern of masses then, clearly, is a complete mystery. Furthermore, it 
has been argued recently that the Jo (1400) should be renormalized down to a very 
wide /o( 1000)0 which would almost certainly imply that / (975) is an exotic. This 
interpretation requires that one identifies the ss member of the nonet which should, 
if GMO is approximately respected in this interpretation, lie somewhere in the 1.8 
GeV range (note that the 0(1700) is a possible candidate). At any rate, it seems 
important to assess whether the 1=1 meson a (980) is a member of the nonet or an 
exotic particle. One of its properties, which was interpreted as an indication for an 
exotic status is the electromagnetic width || (see || for a critical discussion of the 
determination of these widths): 

r(a (980) -> 2 7 ) = (0.24 ± 0.08) / 'Rate(irrj) KeV (2) 

which depends on the unknown irr] rate but, clearly, since this is the only open decay 
channel, one does not expect it to be much smaller than one even though some 
amount of KK decay is certainly observed. In order to appreciate the peculiarity of 
this result, one may compare it with the quark model calculation^]: 

T(a -> 2 7 ) ~ 1.6 KeV (3) 



2 



It is of theoretical interest to understand why the quark model fails by so much in this 
case. We will argue that the amplitude for this decay has to satisfy rather stringent 
chiral symmetry constraints (certainly a surprising fact, at first sight, since no pion 
is involved in that reaction) and this is ultimately why the naive quark model runs 
into trouble. 

3- Chiral correlation functions: 

We intend to use the asymptotic matching technique for the three-point func- 
tion made out of two vector and one scalar current < SVV > which is obviously 
related to the decay amplitude of interest. In order to determine some of the cou- 
pling constants we must discuss also the two-point function < SS — PP >. Let 
us define these objects in a more precise way. The vector, scalar and pseudoscalar 
currents, first, are defined as follows: 

A a A a A a 

jfa) = V>y7^, fs(x) = ^yVs 3 P {x) = #y7 5 ^ (4) 

The couplings of these currents to the resonances which will play a role in what follows 
are then defined (in the chiral limit) as: 

< 0\j;(0)\V b (p, A) >= S ab M v F v e,(X) : < 0\j a s (0)\S b (p) >= 5 ab M s G s 

< 0\3p(U)k b (p) >= S ab B F < 0\f P (0)\P b (p) >= 5 ab M P Gp (5) 

Here, Fq is the pion decay constant and Bq = — < uu > /Fq, which are standard 
notations in xpt . Note that in the pseudo-scalar channel we will take into account 
the pseudo-scalar resonance multiplet (7r'(1300)) in addition to the pion multiplet, 
since it lies in the same mass range as the vector or scalar resonances. Finally, the 
definition of the two and three-point functions is as follows: 

nf P (p) = * J exp(ijtt) < T(j a s (x)j b (0) - j%(x)j b P (0)) > 

VfrM = J d A xd A yexp(tpx + tqy)<Tf a (x)j b i (y)f s (0) > (6) 

Let us begin by considering the two-point function U SP in three successive energy 

ranges: 

a)p — > oo: 

Chiral symmetry implies that the operators which occur in the OPE transform 
in the same way as the operator SS — PP. The important point is that this operator 
transforms non-trivially under the chiral group (like (3, 3) + (3, 3) + ... ). In addition 
there is a discrete symmetry ip L — > — ip L , which implies that the operator of lowest 
dimensionality in the OPE must be of the kind iPl^r^l^r- Ignoring anomalous 
dimensions here, dimensional analysis shows that for large scales: 

hm Uf P (Xp) ~ 1 (7) 

A^oo A^ 
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b)p < 1 - 2 GeV 

In this region, we assume that the lowest lying set of resonances dominate 
the two-point function. In the leading N c approximation the only singularities are 
resonance poles so that, using the definition of the coupling constants fl5|), one obtains 
the following representation: 

Tlf P (p) = 5 ab \^-+ MlGl M|G| 1 



y P 2 p 2 — Mp p 2 — Mg j W 

Now in the very low energy region, 
c)p « 1 GeV: 

Chiral perturbation theory controls the expansion of the correlation functions. 
Up to one loop, one obtains 

II?P(P) = ^ | # " T^L In ^ + Mp) + ...) (9) 



96tt 2 F 2 /i 2 

We notice that, since Fq=0(N c ), the second contribution which is generated by 
the pion loop drops out in the leading N c approximation and, in this limit, Aq([i) 
is a scale independent low-energy parameter. It is not difficult to determine this 
parameter which controls, for instance, the expansion of the pseudo-Goldstone boson 
masses in powers of the quark masses. Up to and including quadratic order, the 
expansion reads 

F 2 M 2 = 2mB F 2 + Am 2 A Q F 2 + ... 

F 2 K M\ = (r + l)mB F 2 + (r + l) 2 m 2 A F 2 + ... (10) 

where 

r = — (11) 

m 

Note that at leading N c order loops are suppressed and consequently, the expansion 
is analytic, xpt also allows one to determine the relation between F and F n to the 
same level of accuracy: 

F 2 2 

7T 

F 2 




Tt = i + — rlii-i) (12) 



from which one can then solve for ttiBq and m 2 Ao in terms of the quark mass ratio 
r and the well known quantities M n , Mk F % , F%. Now the value of r is customarily 
assumed to be r ~ 26, an estimate which may be obtained from (|10D by dropping the 
quadratic quark mass terms. However, an analysis of the violation of the Goldberger- 
Treiman relation in the baryon sector gives instead 6.3 < r < 100]. Here, we will 
assume that r is a free parameter that can vary in the whole range between 6.3 and 
26 and check, eventually, the sensitivity of our results upon that. 



4 



Now, let us perform the matchings. Matching region b) with region c) gives us 
one relation: 

G 2 S - G P = A Fq (13) 
Then matching region b) with region a) gives a second relation: 

M 2 S G 2 S = M 2 P G 2 p + B 2 F 2 (14) 

which is exactly similar to Weinberg's first relation]!]]. Since we can calculate Aq and 
B (or rather mB and m 2 AQ which are RG invariant combinations) as a function 
of r from eqs.(|T0|), we can now determine completely the two coupling constants G$ 
and Gp (note that only the RG invariant combinations mG actually occur in physical 
quantities) which will be needed in the sequel. Before we move to that let us show 
that, already, a non-trivial constraint concerning the scalar meson spectrum emerges 
from the two relations ( |T3"D (|T~]|) : 

Relation (|14D implies that one can parametrize Gs and Gp in terms of an angle 

9 

B F 

M S G S = ^-^ MpGp = B F tan 9 (15) 
cos 9 

Using the first relation we can solve for 9 in the following form: 

M|_ 
tt r 



—% where M n 2 = 

Ml , m 2 Ar 



cos 2 9 = — § where M 2 = — (16) 
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Since cos 2 9 must be positive and smaller than one, M Mp and Ms must satisfy one 
of the two inequalities: 



M P >M S > M Q or M >M S > M P (17) 

Now Mq is a function of r which is easly seen to be increasing and reaches M = 980 
MeV for r = 29.5. This means that for any reasonable value of r the first inequality 
must hold i.e M P > M s . Since M P ~ M n > = 1.2 - 1.3 GeV this inequality is 
compatible with the assumption that Ms = M ao = 980 MeV. On the contrary, it 
would be violated if ao is assumed to be an exotic and the true nonet mass were 
assumed to be 1.4 GeV. 

Let us now turn to the VVS correlation function and, again, go through the 
three energy regimes. 
a)p, q — ► oo: 

It is a simple matter to check that for A — > oo, V Q /3(\p, Xq) scales as 1/A 2 . In 
addition, the leading term of the OPE is the scalar current ipip and, consequently, 
the corresponding Wilson coefficient does not carry any anomalous dimension. This 
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means that the asymptotic behavior is exactly 1/A 2 , and one can as well impose 
matching to the full Lorentz structure of the three-point function V a p. This Lorentz 
structure can be found from the QCD tree graphs to be: 

Vfrip, Q)p^co ~ dabC ^£^ {(P 2 + r2 )(P 2 -I 2 - r 2 )^ + VlW* 
+±p 2 q a qp + 2{p 2 + q 2 - r 2 )p a q p + 2{p 2 + q 2 + r 2 )q a p p ){l + 0{a s )) (18) 

where r = — (p + q). Now in the resonance region: 
b)p, q<l-2 GeV 

We will include the poles corresponding to the lowest lying vector and scalar 
nonets. The VVS Green's function is transverse p a V a p = q®V a p = 0. It can be 
expressed in terms of the two independent tensors with this property: 

P c$ = PpQa ~ P-q Va/3 and = q 2 p a Pf3 + p 2 q a qp - P-qp a qp - p 2 q 2 V a p (19) 

The most general form of the amplitude which matches with the asymptotic behavior 
flT8|) is the following: 



(20) 

One observes that the expression contains a single arbitrary parameter: a. As before, 
it can be determined by matching with the low-energy regime: 
c)p, q^O 

where (forgetting about subleading pion loops) VVS has the following expan- 
sion: 

9c 
5m 

One may recognize here the low-energy parameter c which contributes to the reaction 
77 — > 7r°7r at chiral order larger than four, which was recently studied in [|TTJ ( in this 
reference the parameter is called ofo) and in [[12]. A sum-rule was established for c in 



V$*(p, q) = -—d abc {p p q a - p.qn ap ) + ... (21) 



[12 1 which yields the following estimate 

-(4.6±2.3)1(T 3 

c= 1 

r — 1 



(22) 



The relation of the parameter a to the low-energy constant c is easily determined by 
matching the low-energy region: 

9 

ma = -M^M 2 s c (23) 
o 
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The preceeding discussion has a bearing on the effective lagrangian description of 
the vector and scalar resonances, which is interesting as such, and will provide an 
alternative derivation of the estimate (22) for c. 



4- Effective lagrangian for vector and scalar resonances: 

We will denote the vector and scalar resonaces by and S respectively and 
will assume that they both transform according to the non-linear representation of 
the chiral group (see e.g. ||) as R — > h^Rh'fa). The effective lagrangian can be 
used to compute the three-point function VVS and it has to reproduce the same one 
as found above (p0|) . This uniquely determines the interaction part of the lagrangian 
with two vectors and one scalar fields. Writing C e ff = £s + £v + £vvs, one has: 



C s = (V.SVS - M 2 S S 2 - 2M s G s Ss) 

C v = {--V^ + M 2 V^ + hV^n (24) 



and 



Cvvs ~ F 2 M S G S { M {V " U ~ fvV " v)2S 



-B F 2 [2M V V^S - V^V^ - fyvnS + - fw^) 2 V 2 S\) (25) 

where fy = Fy jMy. (A form which is manifestly covariant is obtained by replacing 
the vector source v^ v by l/2f+ and the scalar source s by l/2x + ). We will use this 
effective lagrangian for two purposes. Firstly, if one is interested in the decay S — > 27 
then there are two terms only which contribute (in this formalism, the coupling of a 
photon to a vector meson vanishes when the photon is on-shell): 

£s 77 = Gt (aSv^ - ^-V 2 Sv, v vn (26) 



This lagrangian differs from the one which was used e.g. in [|TTJ by the presence of 
the second term, proportional to B . This term is needed to enforce the QCD short 
distance behavior even though it is non-leading at low energies. In a sense this is an 
unsusual situation: up to now it seemed possible to achieve compatibility of effective 
lagrangians for resonances with the short distance properties of QCD by retaining 
only terms with a minimal number of derivatives^. Here, we face a counterexample 
to this situation. Both terms contribute to the decay S — > 27 but only the minimal 
term contributes to the low-energy parameter c (or d^). Notice, however, that this 
remark has no practical bearing on the evaluation of the 77 — > 7r°7r° cross-section in 



refjll 



As a second application, the form of C e ff suggests an alternative way to deter- 
mine the parameter a. Let us set the scalar source equal to the quark mass matrix: 
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s = M. = diag(m,m,m s ). The equation of motion allows us to eliminate the scalar 
resonance field: 

S = -^fM (27) 

M s 



Clearly, the lagrangian (pq) then describes propagation of vector mesons in the pres- 
ence of a vector source and the sole effect of the scalar source is to induce ideal-mixing 
and mass-splitting in the vector nonet. One observes that the quark masses act on 
the kinetic energy as well as on the vector mass term. As a consequence, the actual 
meson masses get expressed as a function of the chiral limit mass My and the quark 
masses in the following way 

M, 2 = , - _,ift, ■ ^=5T . < = M; = ■ — -ML 2 ■ — (28) 



l + 2m s a/F v M 2 s + 0{w? q y p " 1 + 2ma/ ' F V M 2 S + 0{m 2 q ) 



and 



2 M V 

Mr * ~ 1 + (m + m s )a/F v M 2 s + 0(m*) (29) 

Note that the terms proportional to i?o in the lagrangian only contribute at order 
O(m^). From ( p8|) one obtains an estimate of the parameter a in terms of the <ft and 
p meson masses. 



n 



ma = MyAfn / v - ( 4r - At) (30) 
y s 2r-l W my V ; 



p 

If one uses the if* and p masses instead one obtains the same result with a 10% 
accuracy. It turns out, in fact, that eq.([30"D reproduces nearly exactly the same 
expression for the constant c (via (E3f)) as the one found in [T3] by saturating the 



sum-rule with the contribution of the vector mesons (in the approximation where 
Fy = F 2 = F 2 = F£ which can be shown to hold, using to first order in the 
quark masses). A remark is in order here. At first order in the quark masses, an 
equally valid way to determine a would be to expand first the denominators in ( P8[) 
and determine a from the mass difference M 2 — M 2 . The result is easily seen to differ 
from the preceeding one by a factor ~ M 2 /M 2 . This essentially confirms the size of 
the error bar in the determination of the constant c in (^2]). 
5- Results and conclusions: 

Now that the unique parameter a is determined it is not difficult, using either 
directly the three-point function (p0|) or the lagrangian form fl26|) , to evaluate the 
amplitude for the decay S — > 27 and to calculate the width. The result is shown in 
the figure as a function of the quark mass ratio r and for three values of the constant 
c, the central one and the two extreme values allowed by the error bar in (p2]) . One 
observes that for the central value the width is in surprizingly good agreement with 
experiment and depends very weakly on the value of the quark mass ratio. Now if c is 
allowed to be larger (in magnitude) then the width is predicted to be also larger and 
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becomes more strongly dependent on r. The width can become as large as 1 KeV if 
r is large, r = 26, and c is at its upper bound. In that case, the width can still be in 
qualitative agreement with expermiment provided the value of the quark mass ratio 
is of the order r ~ 10. 

At this point we should also mention that there is a further uncertainty in the 
calculation, which is associated with the value of Mp ~ M n > (in the figure we have 
taken Mp = 1300 GeV). Clearly, since the 7r' is a rather wide resonance it is not a 
very precise approximation to treat it in the narrow width limit. Despite all these 
uncertainties, we feel that this calculation strongly suggests that there seems to be 
nothing surprising about the small value of the electromagnetic width of the ao, which 
has been reproduced rather naturally. Chiral symmetry was an essential ingredient 
which allowed us to construct the correct form of the three-point function (or the 
equivalent effective lagrangian) and to determine all the parameters. The important 
role of chiral symmetry in this application is at high energy rather than, as usual, at 
low energy. This approach can be extended to discuss the main decay mode of the 
a , ao — > 7r?7 (see fl3| ). In that case, the dependence upon r is rather strong but it 
is not difficult to reproduce the experimental result. All this seem to indicate that 



there could be nothing exotic about the a (980) (as e.g. in the suggestion of refpi 
) which could well be an ordinary nonet member. An interesting question at this 
point is whether the /o(975) is also a nonet member with an exotic mixing angle or a 



completely exotic state. As has been pointed out and discussed in |15[ this question 
could be answered at Da0ne by measuring the decay rates — > a 7 and — > foj- 
These rates can also be evaluated in our approach as a function of the mixing angle, 



the results can be found in 13 . 
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